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It is shown how a chiral Lagrangian framework can be used to derive relationships connecting 
quark-level QCD correlation functions to mesonic-level two-point functions. Crucial ingredients of 
this connection are scale factor matrices relating each distinct quark-level substructure (e.g., quark- 
antiquark, four-quark) to its mesonic counterpart. The scale factors and mixing angles are combined 
into a projection matrix to obtain the physical (hadronic) projection of the QCD correlation func¬ 
tion matrix. Such relationships provide a powerful bridge between chiral Lagrangians and QCD 
sum-rules that are particularly effective in studies of the substructure of light scalar mesons with 
multiple complicated resonance shapes and substantial underlying mixings. The validity of these 
connections is demonstrated for the example of the isotriplet ao(980)-ao(1450) system, resulting in 
an unambiguous determination of the scale factors from the combined inputs of QCD sum-rules 
and chiral Lagrangians. These scale factors lead to a remarkable agreement between the quark 
condensates in QCD and the mesonic vacuum expectation values that induce spontaneous chiral 
symmetry breaking in chiral Lagrangians. This concrete example shows a clear sensitivity to the 
underlying ao-system mixing angle, illustrating the value of this methodology in extensions to more 
complicated mesonic systems. 


Scalars are the Higgs bosons of QCD and induce 
spontaneous chiral symmetry breaking and their inter¬ 
nal structure holds key information about nonperturba- 
tive QCD. However, understanding the light scalar sector 
of QCD has turned out to be considerably more chal¬ 
lenging than other sectors [l]-Q • Particularly the inves¬ 
tigation of the quark substructure of scalar mesons con¬ 
tinues to pose various complications, mostly stemming 
from the underlying mixings among different quark and 
glue components. On the one hand, the experimental 
data on light scalars, which originate from various low- 
energy decays and scatterings, seems to be less certain 
than other better understood channels such as vectors 
or pseudoscalars, while on the other hand, the theoreti¬ 
cal foundation for understanding scalars faces challenges 
within the quark model. The theoretical challenges in¬ 
clude describing the light and inverted scalar mass spec¬ 
trum below 1 GeV, which finds a fundamental explana¬ 
tion in terms of four-quark states described by the MIT 
bag model 3, as well as describing the existing devia¬ 
tions from a quark-antiquark substructure of scalar states 
above 1 GeV, which finds a natural explanation in an un¬ 
derlying mixing among scalars below and above 1 GeV 
[U (other works on mixing include Refs. [il-fl3j|). QCD 
sum-rules also find evidence of an inverted spectrum for 
four-quark states [mUrI . 

Such mixing patterns and global relations are far from 
the focus of chiral perturbation theory [T7I42T] ] which aims 
to probe the near-threshold region in a systematic way. 
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Those formulations of chiral Lagrangians (both linear 
SEMI as well as nonlinear models £5. 91 l22i - [25j ) that 
explicitly include scalar resonances have been shown to 
be effective frameworks for understanding scalars. In this 
approach, the framework is developed on the basis of chi¬ 
ral symmetry, the mechanisms for its breakdown, and the 
model parameters are found from fits to the available ex¬ 
perimental data on scalars below 2 GeV. Particularly, it 
is seen that the underlying mixing patterns among scalar 
states provides an indirect probe into their substructure 
[9, U, 23]. It is shown that the scalars below and above 
1 GeV have considerable mixings, with those below 1 
GeV being predominantly of four-quark nature whereas 
those above 1 GeV being closer to quark-antiquarks. In 
addition, two of them, / o (1500) and /oG710), tend to 
acquire a substantial admixtures of glue p, QCD 

sum-rules have also provided insight into these gluonic 
mixings (see e.g., [271 1421 and review of earlier results in 
Ref. [33j), as well as other methods [dj. 351 ■ 

While the probe of quark substructure of scalars in the 
chiral Lagrangian approach is robust and coherent, it is 
indirect and limited. One of the limitations is the fact 
that the four-quark chiral nonet can be mapped to dif¬ 
ferent quark compositions, each representing different fla¬ 
vor, spin and color configurations (e.g., “molecular” type, 
diquark-antidiquark type). Since these different compo¬ 
sitions transform in exactly the same way under chiral 
transformations, the chiral Lagrangians become oblivi¬ 
ous to differences in internal quark substructure of the 
four-quark chiral nonet. In Refs. EMI , as discussed 
below, it is shown that U(1)a differentiates between two- 
quark and four-quark nonets, but not between different 
types of four-quark substructures. It is therefore desir¬ 
able to supplement these investigations with a method 
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that directly probes the quark substructure. This limita¬ 
tion of chiral Lagrangians can be in principle addressed 
by a linkage to QCD sum-rules. 

Conversely, QCD sum-rules can benefit from chiral La¬ 
grangians. QCD sum-rules exploit quark hadron dual¬ 
ity and connect fundamental QCD to hadronic physics 
I -1(1 . For probing a light, broad and highly overlapping 
state, such as the sigma meson, a simple Breit-Wigner 
shape may not be a realistic parametrization. In ad¬ 
dition, when there are complicated underlying mixings, 
a realistic modeling of the hadronic side in QCD sum- 
rules that includes the effects of such underlying mixing 
phenomena seems necessary. Chiral Lagrangians (such 
as those of Refs. flll4l3j ) where their model parameters 
have been fixed by various low-energy data can provide 
important information on both the resonance shapes and 
mixings. 

The main purpose of this paper is to demonstrate how 
chiral Lagrangians and QCD sum-rules can play comple¬ 
mentary roles in forming a bridge that ultimately con¬ 
nects the relevant low-energy experimental data to fun¬ 
damental QCD. Key parts of this bridge are scale factor 
matrices relating each distinct quark-level substructure 
to its chiral Lagrangian mesonic counterpart. The scale 
factors and mixing angles are combined into a projection 
matrix to obtain the physical (hadronic) projection of the 
QCD correlation function matrix. This bridge is there¬ 
fore mutually beneficial: Chiral Lagrangians provide phe¬ 
nomenological information to guide QCD sum-rules, and 
QCD sum-rules can provide additional methodologies for 
distinguishing between different chiral Lagrangian sce¬ 
narios, deepening our understanding of the chiral La¬ 
grangian approach. 

Although this work is motivated by the scalar chan¬ 
nel, it is extendable to other channels such as the pseu¬ 
doscalar and vector channels. Here, we use the general¬ 
ized linear sigma model of Ref. El to demonstrate how 
connections can be made with QCD sum-rules, but the 
method is more general and can be extended to nonlinear 
chiral Lagrangians such as those in Refs. nm. Also, 
our focus in this work is the relatively less-problematic 
1=1 channel that allows us to demonstrate the connec¬ 
tion without complicating factors. We intend to extend 
the present work to nonlinear chiral Lagrangians as well 
as other channels (1 = 1/2 and 1 = 0 ) in future works. 

We begin by defining our notation. At the mesonic- 
level, we employ the generalized linear sigma model of 
Refs. EM3 which is formulated in terms of two chiral 
nonets M and M' that respectively represent a quark- 
antiquark nonet and a four-quark nonet (a “molecule” 
type and/or a diquark-antidiquark type) underlying sub¬ 
structure. Both chiral nonets transform in the same way 
under chiral transformations but differently under U(1)a: 

M ->U L MUk, M e 2iv M 

M' U L M' U f R , M' —*■ e~ iiv M' . (1) 

The axial charge is the main tool for distinguishing these 
two nonets. Each of these two chiral nonets can be ex¬ 


pressed in terms of a scalar and a pseudoscalar meson 
nonet 


M = S + i(f> 

M' = S' + icj)' (2) 


where the two scalar meson nonets contain the two- and 
four-quark “bare” (unmixed) scalars 


(S\ 4 n+\ (S'\ a'+ 

S = Uo S 2 «° , 5'= I a'~ S'l 

\k k° Si) \ K >- K >° 



( 3 ) 


and similar matrices for <f> and <//. The framework of 
Ref. JTl] provides a detailed analysis of the mixing be¬ 
tween these two “bare” nonets and how that results in 
a description of mass spectrum, decay widths and scat¬ 
tering analysis of scalar as well as pseudoscalar mesons 
below and above 1 GeV. In this picture, the physical 
isotriplet scalars become a linear admixture of two- and 
four-quark components 


/a°(980)\ 

V«o(1450); 


= L~ 


( (Si - Si) IV2 \ 

V(<S'l-Sf) /V2J 


( 4 ) 


where L a 1 is a 2x2 rotation matrix (here we use the 
matrix computed in El)- 

The transformation properties m as well as decompo¬ 
sitions (0 are direct consequences of the assumed under¬ 
lying quark configurations. The two mesonic-level chiral 
nonets M and M' can be mapped to the quark-level chiral 
nonets M^ CD and M lCiCI) . For example, Eq. 0 implies 


( M Q CD )h a ( q R ) b (q L ) a => (S^ 00 )^ oc q a q b . (5) 

To make the connection to the quark world we need 
to make a specific choice for the proportionality factor, 
where here we choose (<S lQCD ) a = q a q b (there is no loss of 
generality or physical consequences associated with this 
choice). The presence of a composite operator in (0 fore¬ 
shadows our connection to QCD sum-rules because they 
are based upon correlation functions of composite oper¬ 
ators [36]. Similarly, M'Q CD can be mapped to quark- 
level composite field configurations. However, in this case 
there are several options, each representing a different 
angular momentum, spin, flavor and color configurations 
for diquark-antiquark combination. Here we do not list 
such quark configurations and only give the specific form 
below that we have used for our example analysis. 

We assume a direct relationship between MQ cd and M 
via a scale matrix Im that adjusts the mass dimensions 

M = J m M qcd . (6) 


Eqs. 0 and 0 imply that under chiral transformations, 
the physical chiral nonet M transforms as follows 

(/m) _1 M -»■ UlVmT'MTJI 


( 7 ) 
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where by multiplying both sides by Im from left 

M I M U L (I M )~ 1 MUl t . (8) 

If M is to transform in the same way as ©, then 

ImUlIm 1 = Ul => [Ul,Im] = 0 . (9) 

Similarly it can be shown (based on chiral tranformation 
of Mt) 


and this is related to the projected physical currents 

J p = L- l I a J QCD (15) 

that probe the physical particles through correlation 
functions 

n^ D 0r) = (0|T [jQ CD (x)jQ CD (0)] |0) . (16) 

Then the projected physical correlators can be defined 


[Ur, Im) = 0 (10) 

which means Im is a multiple of unit matrix, Im = 
x 1, where A (that must be determined) is a scale 
factor connecting mesonic and QCD fields. The quark 
mass factor m q = ( m u + rrid)/2 will ultimately result 
in renormalization-group invariant currents as discussed 
below. Similarly, 


M' = / M 'M ,qcd (11) 


and we can show that Im 1 is a multiple of unit matrix, 
i.e. Im' = y\ 7 T x 1, where A' is the scale factor for the 
M' structures (similar to A, the scale factor A' must also 
be determined). 

At this point the general methodology should be evi¬ 
dent: each distinct substructure (e.g., two-quark versus 
four-quark) requires a separate scale factor connecting 
the meson fields to the QCD currents. Methods for de¬ 
termining these scale factors are developed below, and 
are a central feature of our analysis. To examine the va¬ 
lidity of the connections between chiral Lagrangians and 
QCD sum-rules, we consider a specific example of the 
isotriplets ao(980) and ao(1450), for which 


( a o(980) \ 
V«o(1450); 


= L: 


( (SI - SI) /V2 \ 

/V2j 


L~ l I a ( (S QCU )l-(S QCB )l 

V2 1 v (5' QCD )|-(5' QCD )' 


( 12 ) 


where L a is the rotation matrix for isovectors defined in 
0 and I a is formed out of the scale factors defined for 
the two chiral nonets in (O and (fill) 


f cos 9 a - sin 9 a \ 

\sin 9 a cos 9 a ) ’ “ 


-■a, 

A 3 

0 



(13) 


We can now associate the QCD-level fields and 

S' ^ CD with similar-structure currents that are needed 
for QCD sum-rule techniques, so that the physical cur¬ 
rents can be parameterized in the same way as that of 
states: 


jQCD _ J_ 

V2 


(5Q cd )J- (S^ cd ) 2 2 

QCD^ 1 QCD^j 


( 14 ) 


n£(x) =(0|T [Jf (x)J p (0)] |0) 

= (L- 1 ). l (I a ) lm n%™(x)(I a ) nk (L a ) kj . 

We thus obtain (after a Fourier transform to momentum 
space) the projected physical QCD correlation function 
matrix 

n p (Q 2 ) = f°nQ CD (Q 2 )r°, r a = i a L a (is) 

where T a denotes the transpose of the matrix T a and 
Q 2 = — q 2 . Although (1151) has been developed for the 
specific case of the ao system, its form is quite general 
and can be easily extended to other systems by augment¬ 
ing the matrix I a to include a scale factor for each sub¬ 
structure and expanding the mixing matrix L a accord¬ 
ingly. The key point is that a matrix 7~ a , composed of 
scale factors and mixing angles, can be used to obtain 
the physical projection of the QCD correlation function 
matrix. 

On the other hand, the hadronic contribution to the 
physical correlator is determined via the mesonic fields 
and a QCD continuum 

n fj ( q 2 ) = / d A x e i9 ' x (0|T [Aj(a:)Aj(0)] |0) 

(19) 

+ cont. I . 

The effect of final-state interactions in the nri channel is 
not as significant as this effect in the 7T7 t (nK) channel 
in which the sigma (kappa) is probed, and therefore in 
the first approximation is neglected here. The last term 
represents the QCD continuum contribution inherent in 
QCD sum-rule methods [3jl]. This hadronic correlator is 
related to the projected physical QCD correlator (fTSl) by 
a standard dispersion relation. 

Eqs. mu, © can be used to express some constraints 
on the correlators. In the isotriplet channel case the van¬ 
ishing of the off-diagonal elements of n p (which is a 2 x 2 
matrix in this case) leads to the following sclf-consistency 
condition for II^ 2 CD : 

ttQCD 
li 12 

This expression can be used in various ways. If the off- 
diagonal QCD correlator is known, then this expression 


T a TT 

'll 11 !! 


QCD 


V 


12 


_l_ / T a 17 
i '12 ii 22 


QCD 


T' 

'2 


22 


n~a n~a / ~ra n~'ci 
Ml '99 T '19/99 


( 20 ) 


= Si 


mt; - q* - im a 
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can be used to assess the self-consistency of (or to aug¬ 
ment) the analysis. Alternatively, because off-diagonal 
correlators typically require higher-loop calculations and 
may be unknown (and impractical to calculate), this ex¬ 
pression can be used to eliminate the off-diagonal corre¬ 
lator. This latter approach will be used in our analysis 
of the clq system. 

Following standard QCD sum-rule methodologies, an 
operator (e.g., Borel transformation operator) is ap¬ 
plied to the dispersion relation connecting the QCD and 
hadronic contributions to the projected physical correla¬ 
tors (36). Because the Laplace QCD sum-rules emphasize 
the lowest state, they are not ideally suited to our pur¬ 
poses because they may obscure the mixing phenomena 
central to our considerations. We therefore focus on the 
lowest-weight Gaussian sum-rules that provides similar 
weight to ground and excited states [33, E 


Go(s,r,s 0 ,s t h) = 


\Z4nr 


exp 


4 T 


which is a diagonal matrix 


—ImII H (f)df 

7T 

( 21 ) 




( (G °o )n 


(g|)J 


( 22 ) 


As an illustration of how our proposed approach 
bridges chiral Lagrangian and QCD sum-rule methods we 
now calculate the scale factors A and A' for the isotriplet 
ao scalar meson system. The QCD currents are E, E 


jQCD = M j } Jx = _ dd ) 1^ (23) 

J 2 = (d a ^"frCs T p -a-n- /3) 

6 - (24) 
+ ^S-d^Cj^sp (d^Csp + aO/3 ) - u ^ d 
v 2 

where C is the charge conjugation operator and cot (j> = 
1 /E]. The diagonal terms in (1221) give 

Gfi(s,r) = aAG^° - bBGg CD ( s,t,s «) 


(25) 


G 22 (s,r) = -aBGff D (s,t,s ( 0 2) ^J + bAG%£ D 


( 2 ) 


A = 


COS 2 9 a 


COS 2 9 a — sin 2 9a 


B = 


sin 2 9 a 


_ m q u- 
a a 6 ’ ^ /aaio 


COS 2 9a — sin 2 9a 
1 


(AO 


(26) 

(27) 


where the QCD continuum has been absorbed into the 
QCD Gaussian sum-rules, G^ x and Gf(> respectively rep¬ 
resent ao(980) and ao(1450), and the factor of m 2 is com¬ 
bined with G^ d for renormalization-group purposes. 
Note that each physical sum-rule has its own continuum 
threshold represented by and Sq 2 0 Given an input of 


cos 9 a from chiral Lagrangians El and the physical mass 
and width of the ao states, one can solve (P?5l) for the (con¬ 
stant) scale factors A and A', and develop a procedure 
for the optimized values of the continuum thresholds that 
minimize the s dependence of the scale factors. Expres¬ 
sions for the necessary QCD sum-rules can be found in 
Refs. [H, H, E along with standard values of the QCD 
input parameters (e.g., QCD condensates). We choose 
r = 3 GeV 4 consistent with the central value used in 
Refs. jH,E. In Fig. [Qwe show the s dependence of the 
scale factors for the optimized values of the continuum. 
The remarkable independence of the scale factors on the 
auxiliary sum-rule parameter s clearly demonstrates the 
validity of our proposed methodology. 

The extracted values of the scale factors from Fig. [T] 
A = 115 MeV and A' = 307 MeV, are characteristic of the 
energy scales intuitively expected to emerge in relating 
QCD to chiral Lagrangian approaches. However, A is 
also related to the vacuum expectation value 

(Si) = = 0.056 GeV, (28) 

in remarkable agreement with the chiral Lagrangian pre¬ 
diction of (Si) = 0.061 GeV El, providing strong sup¬ 
porting evidence for our proposed methodology. Simi¬ 
larly, A' is related to the vacuum expectation value 

(S"i) = LSlMM. ~ o.04GeV (29) 

compared with the chiral Lagrangian prediction of 
(S"i) ~ 0.03GeV El- Although this comparison is 
not expected to be as robust because it does not con¬ 
tain renormalization-group invariants, the approximate 
agreement is significant. 


Scale Factor (MeV) 

350 

300 r- 

250 
200 : 

150 
100 1 
50 


s (GeV 2 ) 


FIG. 1. The scale factors A (lower solid curve) and A' (upper 
dashed curve) are shown as a function of s for optimized con¬ 
tinuum thresholds s= 2.2 GeV 2 and = 4.9 GeV 2 . The 
mixing angle cos 9 a = 0.493 of Ref. [llj has been used. 

We have also explored variations in cos9 a from the 
chiral-Lagrangian value cos 9 a = 0.493 El; increasing to 
cos 9 a = 0.6 and decreasing to cos 9 a = 0.4 leads to clear 
increases in the y 2 measure used to optimize the pre¬ 
dicted scale factors. Thus the QCD content of our pro- 
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posed methodology has sufficient sensitivity to discrimi¬ 
nate between different chiral Lagrangian mixing scenar¬ 
ios. As an additional consistency check, we have verified 
that Eq. for the off-diagonal sum-rule G^ 2 CD agrees 
with the estimated order-of-magnitude effects expected 
from the leading-order perturbative result. 

In summary, we have proposed a general method¬ 
ology that connects chiral Lagrangians to QCD sum- 
rules through scale factor matrices relating the two ap¬ 
proaches. These scale factor matrices combined with a 
rotation matrix then provide a matrix used to obtain 
the physical projection of the QCD correlation function 
matrix onto mesonic states. A detailed implementation 
has been provided for the ao(980)-ao(1450) system, and 
the resulting QCD extraction of the scale factors are in¬ 
dependent of the auxiliary sum-rule parameter and are 
in excellent agreement with the quark QCD condensates 
and vacuum expectation values in chiral Lagrangians. In 
this example implementation, we find sufficient sensitiv¬ 


ity to discriminate between different chiral Lagrangian 
mixing angles. Thus our proposed methodology provides 
a powerful bridge connecting the relevant low-energy chi¬ 
ral Lagrangian models and experimental data to the the¬ 
oretical framework of QCD. We expect that this powerful 
synergy between chiral Lagrangians and QCD will per¬ 
mit future progress on more challenging and controversial 
aspects of low-energy hadronic physics. 

ACKNOWLEDGMENTS 

TGS is grateful for the hospitality of AHF, 
Dr. A. Wolfe and SUNY Polytechnic Institute as well 
as partial funding from UUP Individual Development 
Award (2014, 2015) and President Opportunity Fund 
(2014), while this work was initiated and completed. 
TGS and AP are partially supported by research fund¬ 
ing from the Natural Sciences and Engineering Research 
Council of Canada (NSERC). 


[1] K.A. Olive et al. (Particle Data Group), Chin. Phys. C 
38, 090001 (2014). 

[2] S. Weinberg, Phys. Rev. Lett. 110, 261601 (2013). 

[3] E. Klempt and A. Zaitsev, Phys. Rept. 454, 1 (2007). 

[4] R.L. Jaffe, Phys. Rev. D 15, 267 (1977). 

[5] D. Black, A.H. Fariborz and J. Schechter, Phys. Rev. D 
61, 074001 (2000). 

[6] F. Close and N. Tornqvist, J. Phys. G 28, 249 (2002). 

[7] T. Teshima, I. Kitamura and N. Morisita, J. Phys. G 28, 
1391 (2002); ibid 30, 663 (2004). 

[81 M. Napsuciale and S. Rodriguez, Phys. Rev. D 70, 
094043 (2004). 

[9] A.H. Fariborz, Phys. Rev. D 74, 054030 (2006). 

[10] G. ’t Hooft, G. Isidori, L. Maiani, A.D. Polosa and V. 
Riquer, Phys. Lett. B 662, 424 (2008). 

[11] A.H. Fariborz, R. Jora and J. Schechter, Phys. Rev. D 
79, 074014 (2009). 

[12] A.H. Fariborz, R. Jora and J. Schechter, Phys. Rev. D 
77, 094004 (2008). 

[13] A.H. Fariborz, R. Jora and J. Schechter, Phys. Rev. D 
72, 034001 (2005). 

[14] A. 1. Zhang, Phys. Rev. D 61, 114021 (2000). 

[15] T. V. Brito, F. S. Navarra, M. Nielsen and M. E. Bracco, 
Phys. Lett. B 608, 69 (2005). 

[16] H. X. Chen, A. Hosaka and S. L. Zhu, Phys. Rev. D 76, 
094025 (2007); H. X. Chen, A. Hosaka and S. L. Zhu, 
Phys. Lett. B 650, 369 (2007). 

[17] S. Weinberg, Physica A 96, 327 (1979); J. Gasser and H. 
Leutwyler, Annalas Phys. 158, 142 (1984); Nucl. Phys. 
B 250“ 465 (1985). 

[18] G. Colangelo, J. Gasser and H. Leutwyler, Nucl. Phys. B 
603, 125 (2001). 

[19] I. Caprini, G. Colangelo and H. Leutwyler, Phys. Rev. 
Lett. 96, 132001 (2006). 


[20] J.R. Pelaez, PoS CD12, 047 (2013); R. Garcia-Martin, 
R. Kaminski, J.R. Pelaez, J. Ruiz de Elvira Phys. Rev. 
Lett. 107, 072001 (2011); J.R. Pelaez, Phys. Rev. Lett. 
97, 242002 (2006). 

[21] J.A. Oiler, E. Oset and J.R. Pelaez, Phys. Rev. Lett. 80, 
3452 (1998). 

[22] F. Sannino and J. Schechter, Phys. Rev. D 52, 96 (1995); 
M. Harada, F. Sannino and J. Schechter, Phys. Rev. D 
54, 1991 (1996). 

[23] D. Black, A.H. Fariborz, F. Sannino and J. Schechter, 
Phys. Rev. D 59, 074026 (1999). 

[24] D. Black, M. Harada and J. Shechter, Phys. Rev. Lett. 
88, 181603 (2002). 

[25] R. Escribano, P. Masjuan, J.J. Sanz-Cillero, JHEP 1105, 
094 (2011). 

[26] M. Albaladejo and J. A. Oiler, Phys. Rev. Lett. 101, 
252002 (2008). 

[27] L.S. Kisslinger, J. Gardner and C. Vanderstraeten, Phys. 
Lett. B 410, 1 (1997). 

[28] Tao Huang, Hong Ying Jin and Ai-lin Zhang, Phys. Rev. 
D 59, 034026 (1998). 

[29] S. Narison, Phys. Rev. D 73, 114024 (2006) . 

[30] J. Zhang, H. Y. Jin, Z. F. Zhang, T. G. Steele and 
D. H. Lu, Phys. Rev. D 79, 114033 (2009). 

[31] D. Harnett, R. T. Kleiv, K. Moats and T. G. Steele, Nucl. 
Phys. A 850, 110 (2011). 

[32] T.G. Steele, D. Harnett, G. Orlandini, AIP 
Conf. Proc. 688,128 (2004) arXiv:hep-ph/0308074 . 

[33] S. Narison, Nucl. Phys. B 509, 312 (1998). 

[34] P. Minkowski, W. Ochs, Eur. Phys. J. C 9, 283 (1999). 

[35] F.E. Close, A. Kirk, Phys. Lett. B 483, 245 (2000). 

[36] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. 
Phys. B 147, 385, 448, (1979). 

[37] R.A. Bertlmann, G. Launer, E. de Rafael, Nucl. Phys. B 
250, 61 (1985). 

[38] G. Orlandini, T.G. Steele, D. Harnett, Nucl. Phys. A 
686, 261 (2001). 




